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Abstract 

Evolving from turbulent states the 2D fluids and the plasmas reach 
states characterized by a high degree of order, consisting of few vor- 
tices. These asymptotic states represent a small subset in the space of 
functions and are characterised by properties that are difficult to iden- 
tify in a direct approach. The field theoretical approach to the dynam- 
ics and to the asymptotic states of fluids and plasmas in 2D provides 
a considerable extension of the usual perspective. The present works 
discusses a series of consequences of the field theoretical approach, 
when it is applied to particular problems. The discussion is developed 
around known physical problems: the current density profiles in cylin- 
drical plasma, the density pinch in tokamak and the concentration of 
vorticity. 
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1 Introduction 



The fluids and plasma exhibit in two-dimensions a strong tendency to self- 
organisation in the undriven evolution towards stationary states. As shown 
by experiments and numerical simulation, the 2D fluids and plasmas reach 
states of high coherency of the flow, generated by concentration of vorticity 
in few large scale vortical flows. The process is essentially non-dissipative 
since the energy is almost conserved during this process. The presence of 
dissipation is however essential since breaking up of streamlines and recon- 
nection into larger structures is only possible in the presence of irreversible 
resistive-like mechanisms. The vortex merging is the typical process and a 
large number of studies have been done both experimentally and by numer- 
ical simulation. When the initial state is turbulent one may invoke argu- 
ments related with the inverse cascade of energy in 2D but this approach 
has limitted relevance when the process has evolved to the point where the 
number of cuasi-coherent structures is large: the high number of irreducible 
correlations necessary to describe a statistical ensemble of turbulence with 
embedded structures invalidates any perturbative considerations. 

The problem of evolution towards the asymptotic coherent flow states is 
very complex and possibly different approaches must be developed to examine 
different aspects of it. For the states close to the flnal, organised one, the 
fleld theoretical formulation seems adequate. This is confirmed by the purely 
analytic derivation of the sm/i-Poisson equation (describing the asymptotic 
states of the Euler fluid) and by the derivation of a new equation, for plasma 
and planetary atmosphere vortices, with substantial practical conflrmation. 

The particular nature of the physical systems that are investigated is re- 
flected in construction of the Lagrangian of the field theoretical model: for 
point-like vortices we have to adopt an algebraic formulation of the fields, 
namely su (2). For short range interaction we need the Higgs mechanism to 
generate a mass for the photon, and this imposes a particular form for the 
scalar matter self-interaction. Finally, in particular cases it appears possible 
that the su (2) structure is reduced to the Abclian substructure (this is pos- 
sibly similar to the Abelian dominance) and the nature of the extrema of the 
action functional changes significantly. 

In the present work we apply the results of the models developed be- 
fore. Since these models lead to equations describing the stationary states 
(flows) of lowest energy of the physical systems, the application should con- 
sists basically in solving these equations and confronting with experiments 
or observations. However some of these equations are already known and 
eventually they have been derived from different considerations, hke the sta- 
tistical theory. We will make a comparison between our approach and the 
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statistical one. 

The following partial conclusions seem to be supported by the analysis of 
the applications that are presented below. 



1. the field theoretical model of the current profiles in tokamak is com- 
patible with the Liouville equation. Comparison with the model of J. 
B. Taylor gives interesting suggestions for the physical interpretation 
of the FT parameters. 

2. for the Euler fluid we obtain in FT a possible confirmation of the exis- 
tence of a current of vorticity leading to concentration into filaments. 

3. for the 2D plasma in strong magnetic field we obtain patterns of vortical 
flows that confirm previous numerical calculations. 

4. for the LH transtion, we obtain, after a renormalisation of the Larmor 
radius into an effective Larmor radius, profiles of electric fields at the 
edge (in H mode) that are compatible with the experiments 

5. for the density pinch we are able to build a physical picture that is 
consistent with the idea that the pinch of density is due to a pinch of 



6. for the 2D atmosphere we obtain quantitative results that compares 
(very) well with the observations of tropical cyclone. 

7. for the Abelian dominance model the first results show the existence of 
ring-type vortices. 

The various physically relevant cases are included in the following classi- 



vorticity. 



fication 



Current J Abelian 
Euler Non-Abelian 
Superfiuid Abelian 
CHM Non-Abehan 
CHM Abehan 



Alp + exp (^) = 
Aijj + sinh tp — 
Alp = exp {ip) — 1 



conformal invariant no Pg 
conformal invariant no Pg 



n.a. (Minardi?) finite Pg 
non-topological finite p^ 

topological finite Pg 




Alp = {e^ - 1) 
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2 The Liouville equation 



2.1 Applications of the Liouville equation 

There are at least three physical problems that may imply a description in 
terms of the Liouville equation: 

1. the natural current profile in tokamak plasma 

2. the snake of density in JET and other tokamaks 

3. the filament of current density developed on a particular magnetic sur- 
face in tokamak (see Huysmans) 

The last two phenomena seem to result from a concentration of a scalar 
field (density and respectively the current density) and show robustness. The 
last property is an indication that the system has reached an equilibrium that 
is of a purely nonlinear nature, analogous to the relaxed states, therefore they 
could be derived from the extremum of an action functional. 

2.1.1 The current profile 

The statistical theory of J.B. Taylor (and also : Montgomery, etc.) relies 
on the principle of maximum entropy, keeping the total energy and the total 
number of particles constant. It leads to the Liouville equation. 

The constant entropy principle used by Minardi leads to the equation for 
the current and is well verified by statistical analysis of the peaking factors 
for current, density, pressure. 

The statistical studies carried out on a large set of discharges with the 
purpose of testing the prediction of the Turbulent Equipartition theory have 
suggested that the current density is given by the equation Aj + (A^/4) j = 
where A is a constant. Replacing j = Aip and assuming that two space 
integrations are possible without introducing new physical effects, we would 
have the equation 

A-^ + —ip = 
that can be seen as a small ip approximation to 

A^ 

A^/; + — [exp - 1] = 
We just note that this equation is of a similar form as 

Ax = exp (x) - 1 



5 



that describes the superfluid streamfunction in the Abelian-Higgs model. 
Actually these equations are not compatible, except when taking the La- 
grangian multiplier A /4 negative. This suggests that the principle of con- 
stant entropy used by Minardi cannot be derived from the topological theory 
of vortices in superfiuids. 

We will examine below the Taylor approach to the natural current profile 
in tokamak. 

2.1.2 The snake of density 

The snake phenomenon has been observed in JET during the experiments 
of pellet injection and consisted in formation of persistent density perturba- 
tions at rational-g surfaces. These structures persist over several sawtooth 
collapses and are difficult to explain as magnetic perturbations. On the other 
hand there are indications that the tokamak plasma density has an anomalous 
radial pinch, much larger than that of the neoclassical origin. Apparently 
in the class of phenomena one should also include the persistent impurity 
accumulation in laser blow-off injected impurity, observed in experiments in 
TCV. There are several studies of the statistical properties of the correla- 
tions between the peaking factors for density, current density or pressure, 
with plasma parameters and these studies seem to support the idea of tur- 
bulent equipartition of the theromodynamic invariants. However we should 
note that these studies involve quantities expressed as global variables (like 
averages) and they can hide other dependences not immediately obvious. 

We consider the possibility that the particle density behavior (and par- 
ticularly the snake phenomenon) can be connected with the existence of 
attracting solutions of certain nonlinear integrable equations. 

The reason to consider the sm/i-Poisson equation comes from the exist- 
ing proofs that this equation governs the asymptotic states of ideal fluids, 
or, more generally, of 2D systems that can be reduced to the dynamics of 
point-like elements interacting by the potential which is the inverse of the 
Laplacean operator (Jackiw and Pi, Spincanu and Vlad). This equation is 
however obtained when there are two kinds of elements (like positive and 
negative vorticity) and they are of equal number, n+ = n_. Then the sinh- 
Poisson equation is obtained as governing the states with maximum entropy 
of the discrete statistical system at negative temperature (Montgomery et 
ai). Since the equation for the current density mentioned above is derived 
under the assumption of turbulent equipartition, the two descriptions may 
be related. However, the solutions for the unbalanced system of elements, 
a = n^/n_ ^ 1, A(f> = (A^/8) [exp (0) j \fo. — \/aexp(— 0)] have been ob- 
tained numerically (Pointin and Lundgren) and have been shown to have 
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higher entropy and higher stabihty than those of the sm/i-Poisson, precisely 
the characteristics we are seeking for. The hmiting form of the unbalanced 
equation is the Liouville equation, A0 = exp (—0). 



3 Classical physical interpretation of natural 
current profile in tokamak 

The model is simplified in order to exhibit the essential aspect of current 
self-organization in tokamak. The equation is 



In the work of J.B. Taylor (1993) it is shown that a reasonable assumption 
is that the current density should consists of filaments, acted upon by 
as a velocity field and with z as time. The position of a filament is (z) = 
[xi (z) , Hi (z)] and the equations of motion are (all filaments are assumed 
equal jo) 



Vx (J X B) = 



or 



Bo^ + (B^ • V^) J=0 



with 



B_L = -Vip X n 



■ dyi 



1 dH 

Bo dyi 
1 dH 
Bo dxi 



where 



H 



J2JoU (ri,rfe) 



k<i 




In an infinite region, the Green function U (r, r ) of the Laplace operator is 
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The current distribution and the magnetic flux function ijj (streamfunction) 

J = Y,JoHr-ri) (1) 

i 

= jo'^U (r, 

i 

The statistical description of the system. 

The energy of the system is H and it is strictly conserved and a micro- 
canonical ensemble, where the joint probability distribution of the positions 
of filaments is 

p{{r,})^6{E-H{r,}) 
is appropriate. The entropy 



S = —k J dr n (r) In [n (r) 



is a measure of the number of microscopic configurations corresponding to a 
macroscopic configuration n{r). 

Statistical equilibrium is obtained by maximizing S (the entropy) under 
the constraint of energy conservation and fixed total number of filaments 



E (fixed) = jo J dr n (r) ip (r) 
(fixed) = J dr n (r) 



The continuum version of the magnetic flux function (z-component of the 
magnetic potential), ip is obtained from Eq.([T]) 

ip (r) = jo J dr f/ ^r, r j n j 

The equilibrium current distribution obtained by extremizing 

S-/3E--fN 

is 

J (r) = jo(^(r)) 

= K exp [-iSjoip (r)] 

Natural current density profiles. Take ip to be zero on the magnetic 
axis. Then 

V^tjj = Jo exp {—X^) 
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For circular symmetry in a tokamak of radius a, 

2 / 
-0 (r) = - In 1 + a— 
A \ 

where 



Stt 

Introducing the total current /, 

it is found a relation between the peaking factor of the current density, a , 
and the inverse temperature (3 of the current filaments 



1 + q; 

Uniform current (which means that the whole plasma volume is chaotic) is 
obtained for a magnetic temperature = of 

a = or Tjn — >■ oo 

and it is identified a critical magnetic temperature where the totality of 
the current is concentrated into a singular central filament, 

a ^ oo or ^ = 

Values of the magnetic temperatures between and are not accessible. 

Hollow current profiles correspond to negative magnetic temperatures. 
They are only accessible through the infinite value of the magnetic tempera- 
ture, Tm — )■ cxD, which in terms of profiles means that the current passes first 
through a state of uniform distribution. 

To compare with Field Theory we take the two equations 

FT: Alnp = ± — p 

CK 

This is the Liouville equation after the substitution 

p = exp (ijj) 

It has nonsingular and non-negative solutions for p when it has the form 

Alnp + I7I p = 

Then the convention on the choice of signs in the right hand side follows from 
the choice of sign for k: 



9 



1. when K — \k\ > one has to take the — sign in the right hand side, 
such that 

Alnp = j — -p 

c \k.\ 

or ^ 

Alp + (^^^ exp {tp) = 

2. when k — — \k\ < one has to take the + sign in the right hand side, 
such that 

2e2 2e2 
Alnp = H p = —p 

CK c\k\ 

or 

Alp + -r-r exp {^p) = 



^C\K\^ 

Therefore the equation is always the same, whatever is the sign in front 
of the Chern- Simons term. 

In Taylor's theory, we have 

ViV = Joexp(-AV') 



with 



A — 



8na 1 _ .2 



1 + a T, 



JqO? 



m 



We change the variable 

ip ^ ip — —Xip 
and the Taylor's equation becomes 



(-A)VlV'' = Jo exp (V^') 
ViV^' + (^y) exp (V^') = 



or 

'0 o'o 



A Sua c\k\ 
Then we can translate the results obtained by Taylor: 
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1. when the peaking factor a goes to 0, {i.e. the magnetic temperature 
Trn — >■ oo) the current profile is fully relaxed to a uniform distribution; 
This corresponds of vanishing k in the field theory: no Chern-Simons 
is present. 

2. when the peaking factor a goes to oo, [i.e. the magnetic temperature 
reaches the critical value, Tm — )■ T^) the current is strongly concen- 
trated on the axis. This corresponds to infinite value for k in the field 
theory, \n\ — >■ oo: the Chern-Simons term is largely dominating every- 
thing else in the Lagrangian. 

3. Negative magnetic temperature 

are obtained in the Taylor's model when 

q; < 

or, the current profile is hollow. In field theory this corresponds to a 
change of sign of k. But the equation remains the same. The field 
theory starts with a certain sign of k, then the Chern-Simons term is 
suppressed (taking \k\ — )■ 0) (leading to uniform solution for — )■ — oo 
everywhere, while A'^ may remain finite). After that the CS term 
is re-established but with an effect which is opposite to the previous 
regime. 

Everything should be seen as an evolution on the manifold of SELF- 
DUAL states, or solutions of the Liouville equation. The parameter that 
moves the states on this manifold is k. 

The following quantities have dimension of inverse distance squared 

Jo _ _ 1 

A c\n\ 

where p is a distance. This distance will be the natural unit of space-like 
quantities in the problem. For example if our physical problem is localised 
spatially in the disk of radius a, the adimensional space range is 




We note that the space unit p is proportional with k. We can say that the 
passage of the system from a concentrated current profile to a hollow current 
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profile includes a state of strong localisation, where the natural space unit 
is extremely small, which means that different parts of the system are sepa- 
rated and non-interacting (physically this means chaos and uniform current 
everywhere) . 

4 The equation for the velocity of the fluid of 
point-Hke vortices in static configurations 

In the field theoretical models of developed for the 2D current density distri- 
bution and for the vorticity in ring-type structures in fluids or plasmas, it is 
found that the magnetic potential (that carries the interaction between the 
point-like vortices) has spatial components given by the equations: 

A = ± AV X Inp 

where c is a dimensional constant (in the quantum theroies where the objec- 
tive is to describe the Abrikosov Nielsen Olesen vortices, A = he/ (2e)). 

We want to calculate the contribution of the physical velocity in the 
balance equation at stationarity 

(v • V) V = - Vp 

where p is the scalar pressure in the fluid. The physical velocity is just the 
last term in the expression of the magnetic potential A. 
Using the definition of v we have 

(v • V) V = [(AV X Inp) ■ V] (AV x Inp) 

or 

V X In p — )■ 62 X V (In p) 

We extract for the next calculations the physical coefficient A. This will leave 
some apparent incompatibilities in dimensions, but at the end A should be 
re- introduced. 

(v V) V 

- {[e, X V(lnp)]-V}[e, X V(lnp)] 

= [e'''dk{\np)\d,[e^'di{\np)\ 
= e'h^'dk{\np)didi{\np) 
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This is a vector and we have to calculate the two components in plane. For 

j = X, 

[(v-V)vL 

- e^y{e''dk{\np)d^dy{lnp)} 

the summation over / consists of just one term, / = y. We have considered 
the convention = 1. We make the summations over i and k: 

[(v-V)vL 



= e"'dk{\np)didy{\np) 



= e^ydy (In p) d,dy (In p) + e^^d, (In p) dydy (In p) 
= 9.|^[a,(lnp)]^|-9.(lnp)9,^(lnp) 



The other component, j = y 

[(v-V)v]^ 



ey'{e'''dkilnp)dAilnp)} 
{£'^4(lnp)aA(lnp)} 
e^^^^fc (lnp)(9i(9^ (Inp) 



with the convention e^^ = — 1. Now we make the summation over i and k: 

[(v-V)v]^ 

= -e'''dk {In p)did^{]np) 

= -s^^dy (In p) 9,9, (In p) - e^^d, (In p) dyd, (In p) 
= -9,(lnp)a^(lnp) + a,|i[a,(lnp)]^ 

In the expression of the first component we take into account the differ- 
ential equation verified by p 

A In p = ap 

where a is a constant. Then 

dl{^np) + dl{\np) ^ap 

and replace 

(9^(lnp) = ap-a^(lnp) 
13 



in the expression of the x component 

[(v-V)v], 

= d,\^-[dy{\np)f'^-d,{\np)dl{\np) 



d^^[dy{\iip)f 

-d^{\np) [ap-dl{\np)\ 

dx^[dy{\Tip)f 

-d^ {hip) ap 
+d,{\up)dl{\np) 



= 9.|^[9,(lnp)]^| 

- [dxp] 

+a.|^[a,(inp)f 



where 



1 1 
S = -[dy{\iip)f + -[d^{hip)f -ap 

Another expression can be obtained if we substitute 



P 

dip {d.pY 



p p2 



{9xpf ^ {dypf 
P^ P^ 

A(lnp) 

P 

Ap 

ap 

P 
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Then 



1 Ap 3 

— t: TrCtP 

2 p 2 



In conclusion the x component is 



a; 



and the pressure is 



(v ■ V) v] 



3 1 Ap 
^ = 2"^-2T 

3 A /i N 1 Ap 



If we take 
we have 



p = exp 



and 



3 / 2 N 

1 2 

= U V 

2 

where oo and f are normahzed. 

This is the equivalent physical pressure that exists in a fluid for which the 
Liouville equation is fuUfiUed. We used the field theoretical formulation in 
order to calculate the pressure from the equation of momentum conservation, 
at stationarity. 

NOTE 
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We remark that the quantity 

1 2 

looks similar to the square 
where 

In 



as a " vector potential" of a flat curvature space. 



5 The physical pressure for the EULER fluid 

We take the part of the velocity that is determined by the gradient of the 
vorticity density The balance equation at stationarity is 

Po 



and 



Then 



v'^ = Vu = V {-sinhi/j) 
— — cosh-j/'VV' 

(v'^ • V) v"" = cosh i/j ( VV' • V) (cosh V' V^') 
{diip) sinh ip {diip) {djip) + {diip) cosh ip [didjip) 



(v^ • V) v'^ \j — sinh cosh ifj {diip f {djip) + (cosh ip)'^ {diip) (didjip) 
— cosh ip (Vip) ■ [sinh ip ^Vip) (djip) + cosh -0 V (djip)] 

We use the formula 

(v • V) V = V Ty^j - VX (V X v) 
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Then 




V ( — ) - v'^x (V X v'^) 



- Vuj X (V X Vuj) 



[V (sinhV')]' = [coshV' (VV')]' 



or 



for 



Take 



Then 



Po 



-Vp 



Po = P1+P2 
— 2 cosh t/j 



—2 cosh ^0 V 



= (-2 cosh V') -V [(cosh i/jf {Vt/jf] 



This should be expressed as a gradient of a scalar function. 

In order that the expression above to be written as the gradient of a scalar 
function (pressure) we start from thesecond part of the velocity expressed as 



d 



d 
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A typical gradient along this velocity is 

d d 

V • V = V h V 

""dx y dy 

d d d d 

This may be applied in particular on one of the components of this velocity 



d d d d 

d-x ^P^ -P'^d-x^dy ^P' - P'^ dy, 
du! d'^u! du! d'^U! 

+ 



dx dx"^ dy dydx 



d_ 
dx 



1 fdu 



2\dx 



1 /dcu 



2 \dy 



and 



d d d d 

d-x ^P^ -P'^Tx^dy ^P' - P'^ dy, 

du d'^uj ^ du d'^u 

dx dxdy dy dy^ 



d_ 
dy 



1 fdu^ ' 

2 \dx 



1 (du 



2 \dy 



We can conclude that a certain quantity is suggested: 



— ^Vp 

Po 



with the pressure 



p ~ — 



1 (du 



2\dx 



+ 



1 (du 



2 \dy 



i ^p^ - P'^ 



I ^P^ - P'^ 



This pressure is negative, as expected. Since we are approximately at 
self-duality for the Euler fluid we have 



P 



1 fduy ^1 fdu 



2 \dx J 



1 2 

■-(cosh'0) jVV'l 



dy 

dipX^ (dip 
dx J \dy 

2 
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We use the formulas 

(cosh'0)^ = 1 + (sinhV^)^ 

and 
Then 

p ~ — {cosh ijj)'^ I'Vipf 
2 

= -^(l + a;^)v^ 

-'- 2 I 2 2 

= — -V +6^; V 

This part of the pressure of the fluid of point-hke vortices arises from the 
part of the velocity that comes from the gradient of the density of point-like 
vortices. The calculation is also valid in the case of the Charney-Hasegawa- 
Mima (CHM) fluids, i.e. plasma and planetary atmosphere. We see that 
the pressure is most negative there where uj is higher, i.e. at the center of 
the vortex. However, there v is zero. It appears that the pressure is most 
negative somewhere between the maximum v and the center (maximum of 
bS). Since the variation of the a;- velocity v'^ is governed by the force resulting 
from the gradient of the pressure and the pressure has a negative minimum 
somewhere between the radius of maximum azimuthal wind and the center, 
it results that the pressure gradient acts such as to push the "matter" = 
vorticity toward this point. This is the reason for the creation of the cyclone 
eye. 

6 The curvature and the Self-dual state for 
the Euler fluid 

We take two variables identical with those which at SD give us the zero- 
curvature condition. 

where 

A = real constant 
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and calculate (cf. Eq.(59) and foUowings of Dunne) 

K± = d±A:^ - d^A± + [A±, A^] 

From the definitions we have 

0t = ^IE_ + ^IE^ 



Aj^ = aH 
A_ = -a*E 

then 

A+ = - A0 

6.0.3 The expression of the curvature 

Let us take the + sign: 

K^ = 

= d^A--d_A+ + \A^,A-\ 

= a+ [A_ + A0^) - 9_ (A+ - + [A+ - Ac/., A_ + A(/.^] 
= a+ \-a*H + A + 

-(9_ [ai/ - A (01^+ + 02^-)] 

+ [aii- - A + 4>^E_) , -a*H + A + 0^^+)] 

The first two lines are 

-d+a*H + Xd+(t)lE_ + A9+0;E+ 
-(9_ai^ + A5_0iE+ + Xd-(j)^E_ 
- A(a+(/.; + a_0i)E+ 

+A(a+(/.t + 9_02)£^- 
- (a+a* + d-a) H 
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and the last line (the commutator [^_,_,^_]) 

-aa*[H,H] (this is 0) 
+Xa(j)l[H,E_] (this is -2E_) 
+Xa(f)l [H, E+] (this is 2E+) 
+X(l)ia* [E+, H] (this is -2E+) 
-X^Ml [E+, EJ\ (this is H) 
-AVi02 E^] (this is 0) 
+A02a* [£;_, H] (this is 2E_) 
-X^4)2(t)l[E-,EJ\ (this is 0) 
-A^ 020; [£;_,£;+] (this is -H) 

It results 

(2Xa(pl - 2X<j)^a*) E+ 
(-2Aa0t + 2A02a*) E^ 
{-X^,<t>\ + ^ 

[A+,A-]= 2X{a<Pl-a*<P^)E+ 

+2A (-a0* + a*02) ^- 
-A^ (Pi - P2) ^ 

Adding the two equations 

= [A(9+0; + 9_0i) + 2A(a0;-a*0i)]i?+ 
+ [A {d+cPl + a_02) + 2A (-a0* + a*02)] 
(a+a* + a_a) - A^ (pi - P2)] ^ 
This can be written using notations 

K+ = K+H + K+E+ + K+E^ 

with 

K+ = X{d+(f)l + d^(f)^) + 2X{a(f)l-a*(f)i) 
K+ = A(9+0* + 9_02) + 2A(-a0* + a*02) 
K+ = -(9+a* + 9_a)-A'(pi-p2) 
If this curvature should be zero we have 

9+0; + (9_0i + 2a0; - 2aVi = 
+ a_(/)2 - 200^ + 2a*(/)2 = 
- (9+a* + 9_a) - A^ (Pi - P2) = 
This is to be examined as condition for integrability. 
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6.0.4 The expression of the curvature K_ 



Let us take the — sign: 
X_ = 

= d-A+-d+A- + [A-,A+] 
= d. {A+ - A0) - d+ {A_ + A0^) + [A_ + X(j>\ A+ 
= d.[aH - \{(t)^E+ + (t)^E^)] 
-d+ [-a*H + A + 

+ [-a*H + A {(PIE_ + (j);E+) ,aH-X {(P^E+ + 
The first two hnes are 

d^aH - Aa_0i£;+ - \d_(t)^E^ 
+d+a*H - Xd+(f)lE_ - Xd+(t)*2E+ 
= -A(9+0; + 9_0i)£;+ 
-A(9+0t + 9_02)£;_ 
+ {d+a* + d-a) H 

and the last hne (the commutator [^_,^+]) 

-a*a [H, H] (this is 0) 
+X^la[E_,H] (this is 2^_) 
+X(j)la[E+,H] (this is - 2E+) 
+a*A(/)i [H, E+] (this is 2E+) 
-X^(f)l(f)i[E_,E+] (this is -H) 
-AV2</'i [^+,^+] (this is 0) 
+a*A(/)2 [i/, (this is -2E_) 
-AVi02 [^-,^-] (this is 0) 
-AV;02 [E+, E_] (this is //) 

It results 

(-2Aa0; + 2A0ia*)E+ 
(2Aa0* -2A02a*)^- 
(A^i^t - A^0202) -f^ 



-2A (a(/.; -a>i)^+ 
+2A {a0l - a*02) ^- 
+X' {p, - P2) 
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Adding the two equations (jl]) and (jSD 

= l^X{d+4>; + d^^,)-2X{a4>;~a*^,)]E+ 
+ [-A (9+0* + 9_02) + 2A (a0t - a*02)] E- 
+ [{d+a* + d^a) + A' (pi - P2)] H 

This can be written using notations 

K_ = K^H + K{E+ + K^E^ 

with 

= -A (9+0; + - 2A (a0; - a*0i) 

= -A(a+0t + a_02) + 2A(a0t -a*02) 
i^Q- = {d+a* + d^a) + X'{p,-p,) 

If this curvature should be zero we have 

9+0; + 9_0i - 2a0; + 2a*0i = 
a+0* + 9_02 + 2a0* - 2a*02 = 
(9+a* + a_a) + A' (Pi - P2) = 

6.0.5 Powers of products of these curvatures 

The fact that these three expression are NOT zero is the signature that the 
fluid is NOT at self-duality. 

However these expressions only contain the space derivatives and there is 
no possibility to derive a time evolution. 

We also have, for the powers of the curvature, 

tT{El) = 

tr(E!!) = 

tr(M7) = 

tr {HE"^) = 

tr(i7"E!^) = 

trlH^'E":^) = 

tr {ElE"") = tr {E'lE"^) = 1 

tr(ij2") = 2 

tr(/72"+ij = 
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Let us consider the product of the two curvatures 



tr { {K+H + KtE+ + K+E.) {K^H + K{E+ + K^E_) } 

tr {KflKQH^} the trace is 2 

+tr {K;^K{HE+] the trace is 

+tr {K^K^HE.} the trace is 

+tr {K^K^E+H} the trace is 

+tr {KlK:[E+E+] the trace is 

+tr {K^K^E+E^} the trace is 1 

+tr {K^KqE_H} the trace is 

+tT {K^K^ E_ E+ } the trace is 1 

+tr { 7^2^ E_ E_ ] the trace is 



Then 



ir{K+K_} 



or 



iY{K+K_} 

= 2 [- {d+a* + 9_a) - (p^ - p^)] + d.a) + (p^ - p^) 

+ [A (5+0; + + 2A (a0; - a*0i)] [-A (.9+0^ + d^ck^) + 2A (a0t - a 

+ [A + a_02) + 2A (-a0t + a>2)] [-A (5+0; + 5_0i) - 2A (a0; - 

We note that the last two hues are identical 

ii {K+K_} 
= -2[{d^a* + d.a) + \'{p,-p^)f 

-X' [id+<p; + a_0i) + 2 {a<p; - a*<P,)] [{d+<Pl + d.<P^) - 2 {a<Pl - a>2)] 

Now, since we have 

we see that the second paranthesis is the complex conjugate of the first. 

P+0; + 9_0i) + 2(a0;-a>ir 
= (5;02 + d*_4>l) + 2 (a*02 - a0t) 
= 9_02 + 9+01 + 2 (a*02 - «0D 
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Thus this relationship has been verified. We then have 

tr {K+K_} 
= -2[{d^a* + d_a) + X'ip,-p,)Y 

Naturally, we obtain that 

-tr {K+K^} > 

since it is a sum of squares and the equality with zero is precisely the SD 
equations. 

In the powers 

tr {{K-r) 

only the powers of Kq remains and the products (-^^f)" {^2)^- 
Then we have 



tr((i^-)'") = tr{{{d+a* + d^a) + X'{p,-p2))H} 
= 2[{d+a* + d^a)+X'{p,-p2)] 



2n 

2n 



+ ... 



We can however continue this calculation using equations in close prox- 
imity of the Self-duality. 

Prom the first equation at SD 



we derive 



and 



2^-20*01 = 

oz 



dz* 



2a(f)l = 



2^ + 202a* = 

oz 

2^ + 2a0; = 
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with the notations 



from which we have 
and 



d 

d+ = d^ + idy^2— 

d 

d- = d^- idy = 2— 
oz 



d+d_ = A 
d d 



dz* dz 

We can obtain the two potentials at SD 



Then 



We have 



In 01 



„* 

a 




In 02 


a 




ln0; 


2a02 


dz* 


= -9+0; 


2a*02 


dz 


= -5-02 


2a>i 


dz 


--d^<l>. 


2a0t 


- 2^"^^- 
dz* 


-d^ri 



Ki/X = a+0; + a_0i + 2a0; - 2a*0i 

= 5+0; + 5_0i + (-5+0;) - (5_0i) 

= 

X2/A = a+01 + a_02 - 2a0t + 2a*02 

= a+0t + a_02-(a+0D + (-5_02) 

= 
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Xo = -(9+a* + 9_a)-A2 (Pi - P2) 



where we have 



-94 



-9_ln(0O 



9_ 



= -a+a_-[ln(0,) + ln(0t)]-A2(p,-p2) 

which is real. This expression is also zero at SD since it is the equation 

-^Alnp,-A2(pi-P2) = 



-A'0 + A^2sinh'0 = 
2 



The normal choice for A is 



with L the length of the box. 



We NOTE however that it is NOT necessary that the expression of Kq/H 
to go to zero at SD. We just need it to become a numerical factor multiplying 
the local vorticity uj. Because in this case we have 

-^Alnpi - A^ (pi - P2) = pu; 

and this just produces another number in front of u) and a scaling of the 
space axis. 

However it looks better to say that the SD coincides with K = END 



7 A minimizer for Euler 

We can use the expression of the energy, after applying the Bogomolnyi 
procedure, 

E = 2^tr (D_0)) 

which leads to the equation for the states realizing the lowest energy = 
0. 

We have to express the detailed form of the energy. The only factor is 
defined 

L>_ = L>i - iD-2 (6) 
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then 

To continue we express the components of the potential 

A+ = Ax + iAy = aH 
A_ = A^- iAy = -a*H 



Then 



A^ = ]^{a-a*)H 
Ay = ^,{a + a*)H 



Then 



D (f) ^ { - i— i ] E + ( - I I E 

+ ^ («-«*) 02 [^,^-] 

-z^(a + a*)02 [H,E_] 
D-<l> = (^-.^ + 2l(a-a*)0,-2l(a + a*)0,)£;+ 



or 



-^-0 = ^ 2a 01 E+ 

\ ox oy J 

The algebraic ansatz is used for the Hermitean conjugate, 
have 
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where the adjoint is taken for any matrix as the transpose complex conju- 
gated. The change of the order in the commutators is due to the property 
that for any two matrices Ri and R2 the Hermitian conjugate of their com- 
mutator is 

[Ri,R2]^ = (R1R2-R2R1)'' (14) 
= {R2R1 - R1R2) 
— R2R1 — -RJ-R2 

= rIrI 

(* is complex conjugate and ^ is the transpose operators) and we take into 
account that in the expression of (f)' we have already used the Hermitian 
conjugated matrices of E±. 

The Hermitian conjugates of the gauge field matrices are 



4 = l{a*-a)H^^^{a*-a)H 



-l(a* + a)i/t = _J_(a* + a)i/ 



(15) 



Then 



d 



d 



We recall that 
The we have 



a) [,H]--{a* + a) [,H] 



dx ' dy ' 2 
X {<PIE_ + <P*E 



dx 



+ 1- 



dy 



E. 



dx dy 



E^ 



+-{a*-a)4>l [E.,H] 
+ ^(a*-a)0; [E^,H] 
+ [E_,H] 
_l(a* + a)0; [E+,H] 



(16) 
(17) 

(18) 
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or 



1 



The equation becomes 



+- (a* - a) 0t (2E_) 
+^ {a* - a) <P* {-2E^ 
-l(a* + a)cf>li2E.) 

-^(«* + «)0;(-2^+ 



+ 2 



Then 



-{a*-a)(l>*, + {a* + a) 0^ £;+ 



+ 2 



dz^ 



+ 2a(f)2 E+ 



Here we have made use of the identifications 



d_ .d_ 
dx dy 



and 



d ■ d _ 
dx dy dz 



The resulting equations are 

which represent the adjoints of the first set, as expected. 
We can now calculate the energy 

E = — tr f 
2m \ 

-tr 



^4 



2m 
1 

2m 



X 



dx 



dx dy 

Ml . 

-t— 2a 

dy 



,M _ 2a<j)l ] E_ 



(19) 



(20) 



(21) 



(22) 
(23) 

(24) 



dx 



dy 



2a(pl E. 



'1 + 
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dx 



dy 



+ 2a* 02 E. 



2m 



We take 



and we have 



ox ay 



E= — UA 
2m ' 



d 



In 01 - 2a* 



d 
dx- 



01 = 

02 = 



In 01 



The energy becomes 
1 



1 dp,, dx 

+ 1- 



2pi dx- dx- 



+ 



+ 



2 •^S^2 I o *A 

— I— h za 02 



dx dy 

d 



dx. 



In 02 + 2a* 



Pi exp (ix) 



15, . 9x 



2pi dx- dx- 

1 dp2 . dr) 



2a* 



+ P2 



+ 1- 



2p2 dx- dx- 



+ 2a* 



Consider the Sclf-Duahty. 

Now we can take = l/p2 • This is not in contradiction with the fact 
that we try to NOT investigate the self-dual state, but its neighborhood. We 
take 

Pi = — = p = exp (V^) 

P2 



X 



we have 



^=^|exp(^) 



+ _ 2a* 
2 dx- dx.. 



+ exp (--0) 



2 dx^ dx.. 



Or 



or 



^ = ^ [exp + exp (-^)] 



E = - — [exp (i/j) + exp (-■? 
2m 



1^ + _ 2a* 
2 dx.. 

d 



(V^/2 + ix)-2a* 



This form of the energy clearly shows in what consists the approach to 
stationarity and the formation of structure: 
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1. a constant p on the equilines, defined as the zero of the derivatives on 



2. the potentials a and a* become velocities and they become equal with 
the derivatives along the equilines of the angle x- 

The fact that the physical velocity is the derivative of the angle of the 
complex function (p is already known from the Abelian-Higgs model and is 
valid outside the positions of the centres of the vortices. This also shows that 
X must have the same physical dimension like tl), 




and that the phase x is normalized to a quantity that has the dimensions of 



8 Detailed expression of the Euler current 
8.1 Introduction 

We note that in the derivation of the Bogomolnyi form of the energy it was 
not necessary to impose the static states. Then at this moment the states 
may still have a time evolution, assuming that the fluid evolves towards the 
SD 

In this case we can combine the spatial components of the current density 
and inserting the equation written above we get 



"2m 



We return to the expression of the current in the second (gauge-field) 
equation of motion 
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and analogous 



2K{dyAo-doAy + [Ay,Ao\) = -ir 



-2K{d^Ao-doA^ + [A^,Ao]) = -iJ^ 



Now we combine them 

-i{r + iJ^) = 2K{dyAo-doAy + [Ay,Ao\ 
-i (d^Ao-doA^ + [A^,Ao])) 
= 2K{{dy-ida:)Ao 
-do {Ay + iA^) 
+ [Ay-iA^,Ao]) 
2k 

= —{{d:^^idy)Ao 

% 

-do {A^ - lAy) 
+ [A^ + iAy,Ao]) 

J+ = 2K{d+Ao-doA_ + [A+,Ao]) 
= 2K{D+Ao-doA_) 

where we have introduced the notation 

Now we have two expressions for the current density J"*" 



I 

2m 

J+ = 2K{D+Ao-doA_) 



At stationarity 
and 



9oA_ = 

[0, 



4m K 
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This shows that the zero component of the potential of interaction has alge- 
braic content reduced to the Cartan generator 



A, 



■0 



~ H 



The magnitude of Aq is given by the charge p or by the magnetic field B 
(connected via the Gauss law). 

However at stationarity we still have 



since in the expression of A there is a gradient of an angle. 

Then we can say that the system evolves toward this minimum. 

8.1.1 A note on the possible meaning of the equation of motion 
for the LIOUVILLE equation 

The equation for p is very complicated and has the nature of a constraint 
that determines a family of functions p. For some of them the current is 
NOT divergenceless and there is a time variation of p from the equation of 
continuity. But the new function p resulted from advancing in time must be 
again a solution of the second equation = constraint. 

Let us take the expression of the current 

j = -gpA 

where g is a dimensional constant factor. 



A (r) = V X 




P 



We note that 



V ■ A = 



The current is 




The divergence of the current 



IS 



Vj = 



^V- (pA) 

5[(A- V)p + p(V-A)] 
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and using the Coulomb gauge, as shown above, we have 



V-j = -^(A- V)p 
The equation of continuity becomes 

^-5(A-V)p = 



This should look hke 



| + (v.V)p=0 



with V the divergenceless field of velocities associated (but only at self- 
duality) with A. 

We should identify 

v = -gA 

The density p, which for the ABELIAN Jackiw Pi case, or LIOUVILLE 
case, is the 0-th component of the current J'^ , is constant along the trajec- 
tories of the vector field A, which at the self-dual hmit will be v. 

Then p will evolve adiabatically as the lines of fiow of A. If these lines are 
converging to the center then we will have a concentration of matter field. 

8.2 The non-covariant charge of the FT Euler 

The paper 9410065 Dunne identifies the Abelian, non-covariant charges 

go = tr(^^^) 
The component of this charge vector is calculated 

The k component of the charge must be calculated. We take into account 
that 

A, = ^(a-a*)H (25) 
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d_ 

dx 



+ [A\]]{4>,E+ + 4>^E_) 
1 



^E. + l{a-a*)MH,E. 



dx 



- + {a- a*) 01 



E+ + 



dx 



+ (a - a*) 02 



E^ 



For the Hermitean conjugate 



^ + [A\]]{ct>,E^ + ct>,E_) 



dx 



i + (a - a*) 01 



dx 



- + (a* — a) 01 



E_ + 



dx 



+ (a - a*) 02 



E. 



dx 



+ (a* — a) 02 



^4 



The same calculation is made for the y component 

w = + (0i£;+ + 02£;-) 

= ^£;+ + l(a + a*)0i[i/,£;+] + ^£;_ + ^(a + a*)02[i/,£;- 



9y 



i + - (a + a*) 01 



^2 1 



+ T (a + a*) 02 
dy I 



E^ 



^ + -(a + a*)0i E+ + 
dy I 

-p^ T (a + a ) 01 £/_ + 

9y z 



— + T (a + a ) 02 
1 



T (a + a ) 02 



^4 



Using these formulas we can write in detail the charges 



and consider the first term, 
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i 
i 

2m 



— -{a + a ) 01 

dy I 



1 



9a; 



i + (a — a*) 



9a; 



(9a; 



dx 



- + {a- a*) 01 



+ (a - a*) 02 



+ (a - a*) 02 



E_ + 
E_E^ trace is 1 
E+E+ trace is 
E-E^ trace is 
E+E^ trace is 1 



^ T a + a ) 02 



E, 



Then 



dx 



- + [a — a*) 01 



+ 



902 

dx 



+ (a — a*) 02 



I 

2m 



^ 9a; ' ^ 9a; 
Similarly we calculate the second term 



+ ^^2^ + (Pl + P2) (« - «*) 



and obtain 

i 

2m 
i 



tr 



2m 

+ 

+ 
+ 



-tr 



dx 
90t 



i + (a* - a) 0t 



E_ + 



dx 



+ {a* - a) 01 



dx 

(f)iE^E^ trace is 1 



+ (a* - a) 02 



(0i£;+ + 02£;-) 



9a; 
90| 

9a; 



9a; 

which further gives 

^ 2^ 

i 

2m 



+ {a* - a) 01 



+ (a* — a) 02 



+ (a* — a) 02 



')2E-E_ trace is 
b^E^E^ trace is 
(p2E+E^ trace is 1 



90t 
9a; 



+ (a* - a) 01 



01 + 



9x 



+ (a* - a) 



^'^i + ^<^2 + (Pi + P2) (a -a) 
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The preceding results can be combined into 









i 

2m 


^ dx 


+ ax 


+ (Pi + P2) (a - a*) 


i 




+ Ox^^ 


+ (Pi + P2) («* - «) 


i 








~2m 


"^^dx 





and we note that 

0r 



+2 (pi + (o - a*)] 



ax dx 



■01 



a 

= Pi — [2 X phase of 0,] 

ax 
9x 



and 



a02 502 , o 



dr] 
dx 



Then 



= 

m 



, d^ , \( *\ 



We can write for an analogous form 
i 



This charge is ordinarly conserved 

dx^^ 



0. 



8.3 The current of the Euler FT 

The formula for the FT current in the Euler case is (Dunne) 
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We note few differences: the CS in the Euler case is defined with the 
factor 

K 

which we cannot associate to the sound speed since there is no and 
for Eulcr fluid. 

The two similar quantities are 

L = length of the space box 

uJtot = the total amount of vorticity, an invariant for Euler 
1 



L2 



d ru (x, y) 



NOTE 

On the units in the Abelian CS (Abelian Euler) 



^cs 



^ [ X A-K I d'rA^B 



2c I dt 



This will further be integrated over time to give the adimensional action 
functional. 

c T ^ ^ T 
Here units of k, is physically that of speed, then 

[A] = 



The current for /x = /c (space components) is, after the calculations pre- 
sented in more detail for CHM (next Section) 

= + A2 

where 

Afe = -1 {^t (^fc^) _ (^fc^) _ (^fc^t) ^ + ^ (^fc^t) j 

^ -^([*t,[^^,^]] + [^,[^t,^^t]]) 
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8.3.1 The expression of the first part of the current, Ai 

The terms containing space and time derivatives (here the symbol is re- 
placed temporarly by 0) 



where we have to insert 



0t = 0t^- + </';^+ 



This consists of two commutators. 
The first commutator is 



The coefficients of E_E_ and of E^E^ cancel. The result is 

,,(901 ^ , ,,902 



Here we can use 



and obtain 



iY{E+E_) = 1 
tr(i/2) = 2 
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The second commutator in A'l is 



^ — ^- + A — 



90t 



d(i>; 



90t 96* 
-cp^J^E.E^-cj^.^E^E, 
oxk dxk 



^dxk 



E,E_ 



As above, the coefficients of the terms E^E^ and respectively E^E_ cancel. 
The other represent commutators that can be expressed by H: 

[0,9*0*] 



^ Sxjk ^ dxk 
Putting together these results we have 



H 







2 



2 dxk 



(Pi - P2) ^ 



8.3.2 The expression of the second pEirt of the current, A2 

The second line, A2, is calculated in the text xxx. 



We have to give detailed expressions for all components of the current. 
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The X component 

= -'-{-{a-a*){p^ + p^)H+{a*-a){p, + p^)H} 

= 2'-{a-a*){p,+p^)H 
The y component 

K = ~{i{a + a*){p^+p^)H + i{a* + a){p^ + p^)H} 

= 2^(a + a*)(pi + p2)// 

8.3.3 The time component of the Euler current 

This is given by 

= (/.i0*[E+,E_] + 020;[i?_,^+] 



or 



J° = (pi - P2) ^ 



8.3.4 The total expression of the EULER current J'^ 

Finally 



gives 



(9 

2z(a - a*) (pi + P2) - «^ (Pi - P2) 
2(a + oT) (pi + P2) - (pi - P2) 



H 
H 



We note that 



2 ^ ^ 
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8.3.5 Expression of the Euler current at SELF-DUALITY 



At self-duality (and only at self-duality) we can replace the functions a and 
a* that define the potentials A± with expressions of the functions 0^ 2 
0* 2 coming from the first equation of self-duality, = 0. 
We have 



1 
2 

2 2 



2i{a - a*) (pi + p^) 
dij d{2x) 



P2) 



dy dx 
d{2x) 
dy dx 



(Pi + P2) 



l.d 



1 



" = 2 



1 1 

= -2- 

2 2 



1 
2 



We have 



dip 
dx 



-dip ^ d{2xy 

dy 



d 

2{a + a*) (pi + P2) - (Pi - P2) 

(Pi + P2) - (Pi - P2) 
1.9 



dx 
I ^(2x) ' 

ay 



- (pi - P2) 

L[tP + i2tx)]' 2^ 



= (Pi + P2) ^+ [V' - (2ix)] 2 ' 

- (pi - P2) 



= -^^ (Pl + P2)^- 

We NOTE that at SELF-DUALITY we have 

uj — —4 sinh ip 



and it results 



J+ = (Pi + P2) 9+ [ip - {2ix)] - ^id+uj 

d- [ip + (2ix)] - ^id-uj 



J- = -^^(Pi + P2) 

We try to connect this with the pure self-dual state 

Pi + P2 = 2 cosh Ip 

UJ — —4 sinh ip 

d^u! — —4: cosh. Ip (d+ip) 

d-U! — —A cosh. Ip {d-ip) 
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then 



= -i2 cosh — -i2 cosh (2ix) — [—4 cosh ijj ((?+■?/')] 

2 2 2 

= 2 cosh-?/' 

+i cosh tfj (d^ip) + 2i cosh ip {d^ip) 

NOTE 

Assume that 



d+i/j = 
= 

For monopolar vortices with circular symmetry, ip and cu are constant 
along streamhnes 

J+ = iPi+P2)9+X 

J- = (Pi+P2)5-X 
Jo ^ CO 

where 

Pi + P2 = 2 cosh -0 = 2VT+"a^ 

END 



We note that when 00 = the expression + is not zero. 

In the "second quantisation - like" terminology, Pi comes from 0^ and 
represents the creation of vorticity and P2 represents rarefaction of vortices 
then their sum is a sort of total effort in these two actions. They have no 
reason to be individually and shows that in this theory even the vacuum, 
a; = is obtained from an effort of dcnsification followed by an equal effort 
of rarefaction of vortices, such that the final state is unchanged, is void of 
vorticity. 

NOTE 

Let us make a comparison with the structure of the current in the Abelian, 
non-relativistic, CS, 4*'* order (nonlinear gauged Schrodinger eq.) or Abelian 
version for Euler, or Liouville. 

The only source of rotational in the expressions of J^'^ is 



but we also have other terms which are gradients 



In the case of the SD states for Eulcr. we have 



dy dy 
d 



cosh'^ 



= — icosh'?/'^ (i2Y + z/)) 
dy 



The part {p^ + e^Wjip is basicaUy the physical velocity x V-?/' = v. 



The part {pi + P2) — (Pi — P2) looks like a pinch of point- like 



vortices. 
END 



8.3.6 Connection between the Euler FT current and the equations 
of motion of the initial point-like vortex model 

We look for a different form for the current in EULER case. At self-duality 

Pi + P2 — 2 cosh ip 
Pi ~ P2 — 2 sinh ip 



and 



Version 1 



P1 + P2 



2 cosh jjj 



2 cosh -0 



dx 



H 



+ 2— -i— ^ H 
dy dx dx J 

d_ /2^^_ _ ^ 

dx \dx dy J 



H 



We can use the definitions 



d. 



dx + idy 



9_ = d-r — id,, 



2A 

dz* 

dz 



(26) 



Version 2 
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We introduce the vorticity 

Pi- P2 



— 2 sinh 



and we get 



Pi + P2 



1 dp. . 1 d , , 

^ + 2/ - ^^^^ /^i - P2) 

Pi ay ox pi + p2 ox 

_ A (Inp ) + _ . 1 ^(-2^) 

(9i/ ^ 9a; + P2 dx 

'^.di^ ^ d{2x) _ , 1 d{-2uj) 

dy dx Pi + P2 dx 



H 
H 



H 

d{-2uj)' 



d{2x) ^ ^diij\ ^. 1 
dx dy J + P2 



Comment for this NOTE 

This current is composed of two parts 



J 7. 



P1 + P2 

The first term contains the derivative 

5(iV. + 2x) 



dz* 
a+(# + 2x) 



and the first part is practically zero since the phase ip is the streamfunction 
and it is constant on circles. But the phase x should be growing along the 
streamlines. 

The second term is 



d{-2uj) 



Pi + p2 dx 



This looks like the gradient of the vorticity. 
For the Euler fluid we have 



J?) = 



d{-2u) 



—I 



Pi + p2 dx 
1 d 

(2 smhip) 



2 cosh %l) dx 



dx 
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This should be added to the first current 

^ ^ ^ ~\ dx ^ dy ) dx 



9 The 2D plasma in a strong magnetic field 

This is described by the equation 

M + \ sinh {i)) [cosh {ip) - 1] = (27) 

(see however the Second Part for a certain ambiguity in the apphcation of 
the Bogomolnyi procedure, originating from the absence of a topological 
constraint on the residual energy term. This is due to the triviality of the 
first homotopy group of the manifold of the su (2) algebra). 

9.1 The large scale vortices in a 2D plasma with toka- 
mak parameters 

We have obtained large scale vortical flows in the meridional section of a 
tokamak, described by the equation (l27j) . They correspond to previous sim- 
ilar flows identified by analytical methods which are based on the technics 
developed for the Larichev-Reznik modon. 

Numerical solution for L = 307 : monopolar vortex Physical parameters: 



Figure 1: The streamfunction {^p/B) and the velocity, vg{x,y) 

= 0.003 (m), LP'^y^ = a = 1 (m). After normalization L = y = ~ 

330. The unit of velocity is c, = 9.79 x 10 V^e (eV") (m/s) 

Numerical solution for L = 307: dipolar vortex. 

Numerical solution for L = 307: quadrupolar vortex. 

The amplitudes of the flows is in general small (but not incompatible with 
the previous analytic results) and should be corrected for the redefinition of 
the Larmor radius. 



47 




Figure 2: The streamfunction {(fi/B) and the velocity, V0{x,y) 




Figure 3: The streamfunction (f/B) and the velocity, Ve{x^y) 

9.2 Self-organisation of the drift turbulence 

Compare the result from our equation with the result of a statistical/variational 
theory of drift wave turbulence for the Hasegawa-Wakatani equation. 




Figure 4: A theory of self-organization of turbulence (Hasegawa-Wakatani) 
leads to this radial profile for the potential 



9.3 The LH transition 

We obtain as solution of the Eq. fl27|) profiles of velocity with a deep drop 
at the edge. They are similar with the measured profiles observed after the 
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transition L to H. 

We must note that the following results are not exact solutions of the 
Eg. [2l\). They are approximative solutions, or, in other terms, they verify 
Eg. [2l\ ) only when a larger tolerance is accepted in the GIANT code. Never- 
theless, they appear persistently and we can suppose that they correspond not 
to exact minimum action states, but are local minima of the action or they 
may have a slow drift in the space of system's configurations. 
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The soliiion slreamfunclion i4;(x,y) 




The solution streamfunction i)/(x,y) 




Figure 6: Density, setl5_l. p = 1 and L = 411. Streamfunction (f26). 
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The magnitude of the velocity tangent to the streamlines 




Figure 7: density, setl5_l. p = 1 and L = 411. Vg{f27). 



The magnitude of the velocity vector tangent to the streamiines 
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Figure 8: density, setl5_l. p = 1 and L = 411. vq (f28). 
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The solution streamfunction i|/(x,y) 



0,1 ^ 




0.08 ^ 



0.06 



0.04 ^ 



0.02 ^ 



-0.02 ^ 



-0.04 ^ 



-0 06 
500 



500 



-500 



-500 



Figure 9: Density, setl5_5. p = 1 and L = 411. Streamfunction (f29). 

9.3.1 Physical parameters for the H states 

The input parameters are 

a = 1.25 (m) 
Bt = 2.5 (T) 

From this we can calculate the ion cyclotron frequency 



These are fixed parameters. We need to choose a value for the ion tem- 
perature, that will fix the Larmor radius. Take 



eBr 



9.58 X IO^Bt (Gs) 
9.58 X 10^ X 25000 
239.5 X 10^ (s-^) 



T = 2 (KeV) 



then 



c, 



s 



9.79 X lOV^e(e'i^) 
4.38 X 10^ (m/s) 
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The magnitude of the velocity tangent to the streamlines 
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Figure 10: density, setl5_5. p = 1 and L = 411. f^(f30). 



The magnitude of the velocity vector tangent to the streamlines 
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Figure 11: density, setl5_5. p = 1 and L = 411. vg (f31). 
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^ 

= 0.0018 (m) 

This gives 

L = ^ 

Ps 

= 682.35 

The following table shows to what extent we can expect change of the 
effective (adimensional) length of the integration domain 



T (KeV) 


Ps (m) 


L 


0.5 


9.14 X 10-* 


1367 


1 


0.0013 


967 


2 


0.0018 


683 


3 


0.0022 


558 


4 


0.0026 


483 


5 


0.0029 


432 


6 


0.0032 


395 


7 


0.0034 


365 


8 


0.0037 


341 


9 


0.0039 


322 


10 


0.0041 


306 



A case with T = 5000 (eV), where ~ 0.003 and 

L ~ 411 

Prom the numerical simulations, we have the velocity at the bottom of the 
region of negative values, close to the border (from the figure f28) 

ve,iow - 4 X 10"* 

The unit of velocity is 

Cs = 9.79 X 10^^/^l 
~ 7 X 10^ (m/s) 

Then the physical velocity is 

= 280 {m/s) 
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The radial electric field that can generate this rotation is 
or 

Er = 280 X 2.5 = 700 (V/m) 

9.3.2 Calculation of the radial electric field at the edge 

The input parameters are from JET 

a = 1. (m) 
Bt = 2.5 (T) 

Prom this we can calculate the ion cyclotron frequency 



nii 



= 239.5 X 10^ (s-^) 

and 

Cs = 9.79 X lO^y/Te 

Cs 

We must use a single value for Tg instead of a radial profile. Then the value 
of Pg will not be precise. We adopt the point of view that the sonic Larmor 
radius may be modified by various effects leading to 



1 1 1 A Vd\ 



v,= ^ (28) 

We assume that the system evolves to a state where two parameters be- 
come comparable in magnitude: the velocity of plasma rotation u and the 
diamagnetic velocity Vd 

U^Vd 

This is achieved in a structured way, i.e. the various elements behind these 
parameters evolve individually but in a corralated manner to achieve a new 
state of equilibrium. The gradient of the density profile L„ evolves such that 
the diamagnetic velocity Vd increases on a certain region of plasma, which may 
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be connected with the strong gradient at the pedestal. The plasma poloidal 
rotation velocity u increases due to the redistribution of vorticity leading to 
accumulation of vorticity elements toward the center of plasma. It is possible 
that these parameters attain a relative magnitude such that the effective 
Larmor radius p^-^^ becomes sufficiently large. This makes possible a new 
equilibrium of rotation velocity via a new equilibrium vorticity distribution, 
as a quasi-stationary state. We can take several possible values for them 



Vd/u 




L = a/pff 


\'^9bottom 1 




\Er\ {kV/m 


0.2 


1.118 


273 


1 X 10-3 


979 


2.45 


0.4 


1.29 


236 


1 X 10-3 


979 


2.45 


0.85 


2.58 


95 


2.2 X 10-3 


2152 


5.38 


0.95 


4.47 


63 


0.02 


19580 


48.95 



It stops increasing p^^^ for a certain ratio Vd/u. 

See below for a more detailed discussion on this topic. 

9.4 The density pinch 

9.4.1 Introduction 

The vorticity uj = and the particle density n (r, t) are strongly connected 
in 2D tokamak plasma. Then we must look to the intrinsic evolution of the 
vorticity profile since from this we can draw conclusions about the density 
pinch. When we apply the Bogomolnyi procedure to the action functional 
for 2D plasma and atmosphere we find equations describing the stationary 
states which formally would correspond to self- duality, plus an additional 
energy which does not have a topological character. We can assume that the 
equations describe states which are actually quasi-stationary and that the 
additional (or residual) energy term provides a estimation of the possibility 
that the system still continues to evolve. This assumption is valid if we 
can prove that the evolution is slow, which means that the system behaves 
adiabatically. Although this perspective is not clear yet, we will take it as a 
basis to investigate the evolution of the vorticity profiles in tokamak, under 
the field theoretical formulation. 

We first recall that in the model of J.B. Taylor of natural current profiles 
in tokamak it is identified a quantity playing the role of a magnetic temper- 
ature for the statistical ensemble of current filaments. This temperature is 
related with the peaking factor of the profile and it is shown that when the 
magnetic temperature reaches a critical value, T^, the profile of the current 
shows a high degree of concentration, in a form of a single current filament 
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on the magnetic axis. This is vahd for his statistical (maximum entropy) 
theory whose analytical result is the Liouville equation for the flux function. 
However the theory does not provide a physical reason that a possible evolu- 
tion of the system would consist of an evolution of the magnetic temperature 
toward the critical value. 

In the field theoretical formulation of the same physical model (which 
also derive the Liouville equation) the physical parameters are n (the coef- 
ficient of the Chern-Simons term in the Lagrangian) and the electric charge 
e (which is actually the elementary current jo of the discrete model). The 
later can be absorbed into the definition of the potential since the connection 
between the matter field and the interaction (gauge) potential A^^ needs 
not be parametrized in any physical way, being established at the level of the 
discrete model. Therefore the quantity that is still present is k. On the other 
hand, comparing the Taylor model with the field theory model, we find the 
relationship between the parameters, 

\k\ - ^^^^^ ^ 

JqO?c Tjji/T^ — 1 

from which wc see that a very peaked current profile in the form of a con- 
centrated filament, i.e. 

corresponds in field theory to very high values of k. 

While in the statistical model we cannot find an intrinsic reason for 
to approach (from above) the critical value T^, we can look for the presence 
of an equivalent tendency in the field theoretical model, asking when there 
can be reasons for k to increase adiabatically its value. 

In the field theoretical model for the 2D plasma, where k is identified as 
the sound speed 

K = Cg — Ps^ci 

we can assume that k behaves as Pg, since we take the external confining 
magnetic field as constant, S^^* =const. Then an increase of k, can result 
from an increase in p^. 

Physically this is not acceptable when the variation of the temperature is 
not considered, as was one of our basic assumption. Then what can provide 
us with an increase of the Larmor radius? 

A careful consideration of the physical model underlying the field theo- 
retical formulation may clarify this problem. 

The discrete system is based on the idea that elements of vorticity interact 
on distances of the order of the Larmor radius Pg (the inverse of the mass of 
the photon in our theory) and this parameter is taken constant throughout. 



58 



However in the plasma there is another 'field' that has a strong connection 
with the vorticity: the density n (r, t). Due to the compressibility of the ion 
polarization drift velocity, the density is not constant in the CHM plasma 
(as opposed to the Euler fluid case). The spatial variation of the density 
{i.e. Ln 7^ oo) is the cause for the diamagnetic flow with the velocity Vdia — 
PgCg/Ln. The interaction between point-like vortices of the discrete model 
adapts itself to the presence of two velocities: one is the velocity u, simply 
associated with the vorticity, with which the cluster of point-like vortices 
{i.e. the physical vorticity) is in the relationship = V x u. The other is 
the diamagnetic velocity of the plasma, Vdia essentially induced by density 
gradients. The distance of interaction between two elementary vortices is 
modified with the factor 



since this has been proven that affects directly the Larmor radius, replacing 
the physical parameter with an effective one, p^-^-^ , 



It is this quantity that appears systematically when vortical motion is ex- 
amined in tokamak plasma, as illustrated below by the treatment for a drift 
wave vortex. 

Then we should accept that our field theoretical formulation, based on a 
constant Pg can only be valid on small patches, and there, with actually the 
effective Larmor radius. Then we have that the parameter k, the coefficient of 
the Chern-Simons term in Lagrangian, has a "slow" spatial variation coming 
from the spatial variation of p^-^^ , which, in turn is given by the difference 
between the plasma rotation velocity u and the ion diamagnetic velocity Vdia- 
The first conclusion that we can draw at this moment is that an evolu- 
tion of the density profile that leads to an increase of the gradient {i.e. Ln 
smaller) will increase the diamagnetic velocity and then Vdia tends to ap- 
proach from below the plasma rotation velocity u. In consequence the factor 
fl — becomes smaller and the effective Larmor radius increases. This is 



translated in the field theoretical model in an increase of k. If we follow the 
analogy presented above, according to which the increase of k is equivalent 
to Tjn \T^, we find that there is an enhanced clusterization of the elements 
of vorticity, with a possible evolution toward a single filament in the center. 

This process is not stationary since it has a positive feedback: the clus- 
terization of the vorticity toward the center drags the density (basically from 
Ertel's theorem, but the process is more complex) and the evolution of the 





u 
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density increases the gradient and consequently the diamagnetic velocity. 
The effective Larmor radius increases and the parameter k increases which 
still enhances the clusterization process for the vortical elements. 

However there is a limit to this process, coming exactly from the same 
reason for which the vorticity and the density evolve together and exhibit an 
inward pinch in tokamak. 

The limit consists of the fact that the factor 1 — Vdial u ends up by sat- 
urating the process of pinching, due to a too large increase of the effective 
Larmor radius. 

When the effective Larmor radius is too big, the interaction between the 
elements of vorticity is no more of short range but can be considered of long 
range. Or, this is the case of the Euler fluid, where the range of interaction 
is Coulombian {i.e. In). For the Euler fluid there is no compressibility of the 
background density, the density and the vorticity are decoupled and the den- 
sity cannot follow the vorticity. The compressibility of the ion polarization 
drift is proportional with the inverse of the square of the effective Larmor ra- 
dius and this diminishes accordingly. The fluid becomes less Hasegawa-Mima 
and more Euler. Actually the state where Vdia = u corresponds exactly to 
the two fluids: Hasegawa-Mima (because the density is adiabatic and then 
Vdia = u) and Euler (since there the range of interaction is infinite, Coulom- 
bian). But the directions from which the two fluids arrive at this state are 
completely different. 

This intuitive representation of the density pinch (via the vorticity pinch) 
may be useful. However the field theory introduces an additional factor in 
this model: the fact that intrinsically the quasi- self -dual solutions have an 
adiabatic evolution related with the existence of a residual energy which is 
not of a topological nature. 

For the Eq. fl27p the residual energy puts severe limitations to the increase 
of the normalized value of the streamfunction. Any ip that is greater than 
~ 3 is strongly inhibited by a severe increase of this residual energy. Or, 
the concentration in the center of the vorticity is in general accompanied by 
an increase of ip. We should note however that for a normalized domain of 
integration L = L^^y^ / p^J^ of the order greater than ~ 10 the value of ip is 
already smaller than 1, and the limitations imposed by the residual energy 
are far. Moreover, the graph of this residual energy shows that there is a 
favorable slow evolution for the small values of ip since close to ip = Q the 
energy is decreasing with increasing ip. 

We still have to work on that. 
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9.4.2 The role of the effective Larmor radius p^J^ 

In Stationary vortices drift waves (Nycander) it is derived the equation 



where 

and the normalizations are usual: Pg , VLi and e/T^ for space, time, potential. 
Prom here we derive, taking constant temperature 

In'n 1 11 dnn 

+ = , , , + 1 



uriQ T {x) yphys ^£±j no drP'^y^/p^ 



Qi \nodr J u 



= 1-^ 



where we defined 

such that when 
we have 



CsPs 1 



Ln<0 



Vd>0 



The two velocities, and w are in the same direction (positive) 
when the density decreases radially, which is normal in tokamak. 

Therefore we have that 

l-^>0 

u 

if the velocity of plasma rotation is higher than the diamagnetic rotation. 
When the diamagnetic rotation becomes comparable with the plasma rota- 
tion velocity, just slightly lower, this factor is almost zero and the effective 
Lamor radius becomes infinite. 
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10 Conclusions 



This work is intended to provide arguments in favor of validity and usefulness 
of the Field Theoretical approach to the description of fluids and plasmas in 
the evolution towards asymptotic, coherent structures. 

The basis material consisting of definition of the Lagrangians, equations of 
motion and Self-Duality states are presented elsewhere and is not mentioned 
here. The applications are only shortly presented. 

We focused on various instruments, like the field-theoretical currents, 
energy and action functional in close proximity of the self-duality. We have 
noted that there are important aspects that can be examined using these 
instruments, specific to the field theory: 

• the existence of the current of point-like vortices, which can explain the 
concentration of vorticity 

• the slowing down of the evolution in the cuasi-asymptotic state, where 
the number of vortices is still not the minimum that is obtained as the 
pure extremum of the action at self-duality 

• the importance of parameters like the coefficient of the Chern-Simons 
factor, in the adiabatic change of the regimes; 

• the role of the effective spatial dimension, which is adjusted sponta- 
neously via the gradients of the vorticity field. 

Much remains to be investigated but there are already reasons to develop 
the field theoretical formulations in a large class of fiuid and plasma problems. 
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